Persistent currents in diffusive metallic cavities: Large values and anomalous scaling 

with disorder 
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The effect of disorder on confined metalhc cavities with an Aharonov-Bohm flux fine is addressed. 
We find that, even deep in the diffusive regime, large values of persistent currents may arise for a 
wide variety of geometries. We present numerical results supporting an anomalous scaling law of 
the average typical current < Ityp > with the strength of disorder w, < hyp >~ with 7 < 2. 
This is contrasted with previously reported results obtained for cylindrical samples where a scaling 
< hyp >~ has been found. Possible links to, up to date, unexplained experimental data are 
finally discussed. 



Modern advances in nanotechnology made possible to 
manufacture systems in the submicron range. For those 
structures and at temperatures of few mK , the quan- 
tum mechanical coherence of the electronic wave func- 
tions dominates the transport of charge [Q . 

Recently, billiards with rough boundaries have been 
considered in order to analyze the effect of generic bound- 
ary deformations on the spectral and eigenfunction prop- 
erties of ballistic samples |^^. In addition, during the 
last years, the magnetic response and the persistent cur- 
rents of integrable and chaotic billiards in the presence of 
a magnetic flux have been the object of many research pa- 
pers [|Ml • Most of them were inspired in the experiments 
of Ref. in which the magnetic susceptibility and the 
persistent currents of mesoscopic semiconductor samples 
were measured. Continuous billiard systems posses an 
energy dispersion relation quadratic in the wave number 
k and are reasonable models to describe semiconductor 
quantum dots in the ballistic regime. The situation is 
qualitative different for metallic systems. The breaking 
of the continuous translational symmetry by the lattice 
leads to a band structure and to a dispersion relation in 
k which is equivalent to that of the free electron model 
only in the large wavelength limit. Moreover, the metallic 
regime is characterized by a half filled band and a typical 
electronic wavelength of the order of the lattice constant. 
Also and unlike high mobility heterostructures, a realis- 
tic model for a metallic sample should contemplate the 
effect of the elastic scattering with random impurities. 

The observed persistent currents (PC) in metallic sam- 
ples 10| have been typically much larger than the pre- 
dicted theoretical estimates that take into account the 
effect of static impurities and much effort is still devoted 
to understand the nature of the discrepancy between the- 
ory and experiments [Tl| ] . In a recent work we have shown 
that by effect of confinement, large values of PC can arise 
in the half filling regime for clean metallic cavities of dif- 
ferent shapes 

In this letter we address the effect of bulk disorder in 



confined metallic systems threaded by a magnetic flux. 
We are mainly concerned with the interplay between con- 
finement and disorder effects and in particular with the 
sensitivity of the persistent current to the strength of 
disorder w. Our main result is that even deep in the dif- 
fusive regime large values of persistent currents can still 
be obtained for a wide variety of sample shapes. 

We employ the spinless tight-binding hamiltonian de- 
fined on a square lattice and introduced in Ref. p^ ] , with 
the additional inclusion of diagonal disorder. Once the 
profile of the sample is defined, the on-site energies are 
taken randomly distributed in the interval [— w/2,w/2]. 
We consider non-single connected loops of various shapes 
with up to 3 000 sites and values of w ranging from 0.02 
to 5 in order to study the transition between the ballis- 
tic and the diffusive transport regimes. In the following 
all the lengths will be in units of the lattice parameter 
a, and the magnetic flux a will be measured in units 
of ipo = hc/e. An isolated non-single connected sample 
threaded by a magnetic flux a carries a persistent cur- 
rent /, even in the case of multiple scattering by static 
impurities. For a sample containing N non-interacting 
electrons, the total current can be calculated at zero tem- 
perature as / ~ —dE/da with E the total energy of the 
system In the absence of disorder, the typical cur- 
rent hyp = \J /2 da gives a measure of the magnitude 
of the current and is a (fluctuating) function of TV. We 
denote its value at half filling /q. 

For a finite value of w, we define the average typi- 
cal current < hyp >= \J < Ityp >Sn performing averages 
in a small energy window containing 1 << Sn << N 
levels. It was previously checked that this procedure is 
equivalent to average at the same time on the disorder 
configuration (w) and on the number of levels (Sn) that 
is, < ... >^^Sn=< ■■■ >Sn 

Before considering disorder effects, we summarize the 
main results obtained for clean geometries. In Ref. |l^] 
we show that for annular and square loops the corre- 
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spending density of states (DOS) has a local peak close 
to the band center {E = 0) at zero magnetic flux. In the 
absence of surface roughness (SR), the peak consists of 
border states that have in general large values of mean 
angular momentum. These states contribute to enhance 
in one or two orders of magnitude the value of lo with 
respect to values obtained for geometries with the same 
number of sites but without border states. The peak in 
the DOS is robust to changes in the size of the sample 
while preserving its shape. In addition, we here present 
numerical evidence supporting the conjecture that the 
peak in the DOS is present when the symmetry of the 
underlying lattice conforms a symmetry subgroup of the 
confining potential. This is illustrated in Fig. |^ , where 
the spectra for four different samples are shown close to 
half filling: (a) square SL, (b) annular AL, (c) octagonal 
OL and (d) pentagonal PL loops jl^. In Figs. 0(a), (b) 
and (c), a bunch of degenerate levels is observed at zero 
flux for E = 0. These levels evolve paramagnetically for 
finite a, generating large values of PC . Due to absence 
of crossings between occupied and empty levels close to 
half filling, Ityp increases monotonically with the num- 
ber of particles [|l2|. The PL, Fig.^d), corresponds to a 
generic geometry where the mentioned symmetry condi- 
tion is not fulfilled. There is no bunch of states at zero 
flux, and therefore no enhancement of the PC near half 
fiUing is expected. 

When the confining potential has not the square or 
rectangular rotational symmetries our samples posses 
surface roughness (SR). This is the case for the AL and 
the OL in which, besides the border states that carry 
large currents, we find localized states in the bunch of 
degenerate levels at zero flux. These are states trapped 
by the SR and are not sensitive to the external mag- 
netic flux. The flat lines on the top of the spectra of 
Figs. |l| (b) and (c) correspond to these states. A re- 
markable feature is that the charge distribution for the 
border states seems to be highly concentrated on classi- 
cal periodic trajectories that bounce off the borders of 
the loops. We find, as in quantum billiard models, bor- 
der states that look like resonant states associated to a 
continuous set of non isolated periodic orbits In 
the left panel of Fig.^ (a) we show for the square loop 
SL, the contour plot of the absolute squared value of the 
wave functions (|^(a;,y)p) for the 20*'' state below half 
filling for a — 0.025. In the right panel of Fig. || (a) 
the family of classical orbits that supports the quantum 
state is drawn for comparison. Figj^ (b) shows for the 
annular loop AL the contour plot for a quantum state 
with energy close to zero (that is belonging to the peak 
in the DOS) and with the charge distribution extended 
along a diamond-like orbit that goes around the sample 
with large mean angular momentum. As a consequence 
of the SR there is not a complete family of trajectories 
that supports the quantum state like in the square loop. 

In summary we have found, as a consequence of the 



symmetry properties of the confining potential and the 
underlying lattice, clean metallic geometries with large 
values of PC at half filling. We name these kind of sam- 
ples symmetric loops. However, a realistic model for a 
metallic cavity should contemplate the inclusion of bulk 
disorder in order to mimic the diffusive transport regime. 

Sustitutional disorder introduces in the model a new 
energy scale that competes against the existence of bor- 
der states carrying anomalously large currents. In order 
to study disorder effects on the PC it is necessary to 
determine the range of values of disorder w that corre- 
spond to the diffusive regime in our samples. We begin 
analyzing the scaling of < Ityp > with w for the SL. 
According to elementary scattering theory a scaling as 

< Ityp >^ w~'^ is predicted in the diffusive regime Q. 
In Fig. ||(a) we show with square symbols the Log-Log 
plot of < Ityp > at half filling as a function of w for 
the SL in which the total number of sites is = 2688. 
We clearly observe a range of values of w between 1 
and 4, before the exponential damping characterizing 
the onset of the localized regime, in which the scaling 
is < Ityp >~ 'w~'^ , with 7 = 1.98 obtained from the 
best fit. Therefore, this range of w should correspond to 
the diffusive regime for other geometries containing the 
same number of total sites. Indeed we have obtained in 
the same range of values oi w a similar scaling law for a 
rectangular loop with N = 2424 (both rectangles in the 
loop having the same aspect ratio), even though this loop 
is not a symmetric loop (in the sense defined above). Nu- 
merical calculations performed on cylindrical geometries 
satisfied the scaling law < Ityp >^ w^^ p^ , p^ . 

On the other hand for samples with SR like the annu- 
lar loop AL, the dependence of < Ityp > at half filling 
with w shows a non monotonic behavior that changes 
the overall value of the exponent to 7 7^ 2 (see the cir- 
cles in Fig. 11(a)). We argue that by effect of the bulk 
disorder (and before the onset of the localized regime) 
the states trapped by the SR are unpined being able to 
carry a finite PC. Moreover, there is no reason a priori 
to expect the same degree of response to a perturbation 
for all the trapped states. This explains the step-like 
behavior of < Ityp > with w observed in Fig. ^(a) for 
the AL. This non-monotonic response of the PC to dis- 
order translates, in the case of the annular sample, to 
an average exponent 7 « 1.3, significantly lower than 
in the case of the square loop SL. As a consequence, 
as can be observed from Fig. ^(a) there is range of val- 
ues w where the PC in diffusive annular loops are larger 
than those in square loops, despite the SR present in 
the former samples. This overall change in the exponent 
characterizing the diffusive regime does not seem to be a 
particular attribute of samples with SR. For other non- 
integrable geometries without SR, like rectangular loops 
where the inner an outer border have different shapes 
(chaotic geometry), we also find an anomalous scaling 

< Ityp >~ w~'^ with 7 < 2. Nevertheless in the absence 
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of SR we do not observe the non-monotonic response 
characteristic of samples with SR. Fig. ||(b) shows the 
Log-Log plot of < Ityp > as a function of w for a non- 
symmetric rectangular loop formed by an inner rectangle 
of sides li — 13 and I2 — 8 and an external rectangle of 
sides Li = 53 and L2 — 48, respectively. The obtained 
best fit is < Ityp >~ w~^-^^. We stress that in the ab- 
sence of disorder this geometry does not present border 
states because the symmetry condition that leads to the 
peak in the DOS is not satisfied. In other words, Iq for 
a clean non-symmetric loop (see Fig. ^ dotted line) is 
much smaller than Ig for clean symmetric loops (see in- 
set of Fig. I). Even more, Ityp for diffusive symmetric 
loops can be of the same order of magnitude than Ityn 
for a clean non-symmetric loop, as it is shown in Fig. |j 
for a clean pentagonal loop PL and an annular loop AL 
with u; = 1. 

In order to determine for a given diffusive sample the 
value of < Ityp > at half filling, two factors must be 
considered: the value of the current in the clean system 
(/o) and the correction factor that takes into account the 
effect of disorder. For generic non-symmetric systems, 
the exponent 7 that takes into account the diffusive cor- 
rections is 7 < 2 but, as we previously explained, no 
enhancement of the value of Iq is expected. On the other 
hand, in clean symmetric loops without SR due to the 
presence of border states, large values of Iq can be ob- 
tained, although in the diffusive regime the exponent is 
7 = 2. It is for clean symmetric loops with SR, that 
we find an optimal interplay between geometrical effects 
and bulk disorder. That is large values of Iq compared 
to those of generic samples, and an anomalous scaling 
< Ityp >^ with 7 < 2. Therefore these systems 
are candidates to have the largest values of persistent 
currents even deep in the diffusive regime. Also, we ob- 
tain that in symmetric loops (with or without SR) the 
response of the system to the external fiux is always para- 
magnetic, even in the difussive regime. 

The novel effects described in this work can be related 
to up to date unexplained experimental data on PC in 
mesoscopic metallic rings. In Ref. [|j rectangular and 
circular Au loops with a large number of channels were 
designed and the PC measured. The corresponding sig- 
nals were paramagnetic and more than an order of mag- 
nitude larger than predicted values that include the dif- 
fusive correction factor. In addition, in that experiment 
the reported value of the signal for the annular loops was 
5 times larger than for the rectangular loop (although 
the estimated theoretical values were roughly the same). 
We believe that in order to model the samples employed 
in this experiment, geometries with decoupled degrees of 
freedom (like the cylinders employed in many theoretical 
papers) are not appropriate. We present numerical ev- 
idence showing that by breaking the translational sym- 
metry in all directions, it is possible to create at finite 
magnetic flux border states near half filling. This states 



coherently contribute with large values of PC. As we have 
previously shown, the confinement in general gives rise to 
geometrical disorder (SR) that competes with the susti- 
tutional disorder. This results in a non-monotonic behav- 
ior of < Ityp > with w and an overall exponent 7 < 2. 
Therefore, we can obtain for these systems larger values 
of PC than in previous models. Moreover, we find that 
the symmetric annular loops can have in the diffusive 
regime a PC as larger as in generic clean systems. In 
addition we show that in generic non-integrable geome- 
tries without SR the exponent 7 is also < 2, although 
the response of < Ityp > with w is monotonic due to 
the absence of trapped states in the clean regime. This 
anomalous scaling of < Ityp > with w could be a quantum 
signature of the classical non-integrability of the loops. 
The absence of symmetries in the clean system makes 
the effect of the sustitutional disorder not such strong as 
in the case of systems that in the clean regime preserve 
some rotational symmetry. Indeed, for clean integrable 
geometries, like the squares and symmetric rectangular 
loops presented in this work, and also in cylinders 
the exponent is 7 ~ 2. We believe that due to recent 
technological developments it should be feasible to de- 
sign samples with specific symmetries in which our pre- 
dictions could be confirmed. 
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FIG. 1. Energy levels close to the band center as a function 
of a for: (a) SL in which the external (internal) square has 
a side Li — 53(L2 — 13), (b) an AL loop with external (in- 
ternal) radius R = 30(r = 3); (c) an OL loop with the inner 
(outer) octagon inscribed in a circumference of r = 3{R = 32) 
and (d) a PL with the inner (outer) pentagon inscribed in a 
circumference of radius r = 3{R = 35 ). 



FIG. 2. (a) Contour plot of the absolute value square of the 
wave function for a border state of the SL for a = 0.025 (left 
panel). Sketch of the family of classical periodic orbits that 
support the quantum state (right panel), (b) diamond-like 

border state for the AL (left panel) toghctcr with the as- 
sociated classical orbit (right panel). See the text for more 
details. 



4 



A 

p. 

bD 
O 



-2 

A 

p- 

O 




(b) 





Log w 



FIG. 3. (a): Log-Log plot of < Ityp > as a function of w 
for the AL (circles) and for the SL (squares). The dashed line 
shows the best fit < Ityp >~ m"^ *'* for the SL. (b): Idem but 
for a non-symmetric (chaotic) rectangular loop (diamonds) 
with the dashed line showing the best fit < Ityp >~ w~^'^^. 
See the text for details. 
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FIG. 4. It yp as a function of the rescaled number of levels 
A''* = n/N for the annular loop AL in the diffusive regime 
1^ = 1 (thick solid line) and for the clean pentagon PL (thin 
dotted line). Both curves are almost superimposed. The inset 
shows for comparison Ityp as a function of N* for the clean 
AL. 
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